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A one dimensional non-equilibrium stochastic model is proposed where each site of the lattice 
is occupied by a particle, which may be of type A or B. The time evolution of the model occurs 
through three processes: autocatalytic generation of A and B particles and spontaneous conversion 
A ^ B. The two-parameter phase diagram of the model is obtained in one- and two-site mean 
field approximations, as well as through numerical simulations and exact solution of finite systems 
extrapolated to the thermodynamic limit. A continuous line of transitions between an active and 
an absorbing phase is found. This critical line starts at a point where the model is equivalent to 
the contact process and ends at a point which corresponds to the voter model, where two absorbing 
states coexist. Thus, the critical line ends at a point where the transition is discontinuous. Estimates 
of critical exponents are obtained through the simulations and finite-size-scaling extrapolations, and 
the crossover between universality classes as the voter model transition is approached is studied. 
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I. INTRODUCTION 

The phase transitions exhibited by stochastic models 
with absorbing states have attracted much attention in 
recent years, particularly in order to identify and un- 
derstand the aspects which determine the universality 
classes in those models. Most of these models have not 
been solved exactly, but a variety of approximations al- 
low quite conclusive results regarding their critical prop- 
erties. Stochastic models are, of course, well fitted for 
simulations, but closed form approximations and other 
analytical approaches have also been useful in investigat- 
ing their behavior Q. 

One of the simplest and most studied model of this 
type is the contact process (CP), which was conceived 
as a simple model for the spreading of an epidemic and 
proven to display a continuous transition between the 
absorbing and an active state, even in one dimension 0- 
Actually, it was found that the CP is equivalent to other 
models such as Schlogl's lattice model for autocatalytic 
chemical reactions |3] and Reggeon Field Theory (RFT) 
The CP belongs to the direct percolation (DP) uni- 
versality class, together with others models such as the 
Ziff-Gulari-Barshad model of catalysis 01 and branching 
and annihilating walks with an odd offspring The 
DP conjecture states that all phase transitions between 
an active and an absorbing state in models with a scalar 
order parameter, short range interactions and no conser- 
vation laws belong to this class This conjecture was 
verified in all cases studied so far |g|. 



Here we study a generalization of the CP, with an ad- 
ditional parameter, so that the CP transition point be- 
comes a critical line. This model is similar to the model 
proposed by Hinrichsen j§| in the particular case where 
his parameter q is set equal to one. Since the symme- 
try properties of this generalized model are the same of 
the CP, it is expected that this critical line should be- 
long to the DP universality class. However, at one point 
of this line the model is equivalent to the zero tempera- 
ture Glauber model [10], also called the voter model fvi\ . 
which displays a spin inversion (or particle-hole) symme- 
try and therefore belongs to another universality class 
(the compact directed percolation (CDP) class). Thus 
the critical line in the phase diagram of the generalized 
model starts at the CP model and ends at the voter 
model, a crossover between the two universality classes 
being observed. 

In section m we define the model and show its equiva- 
lence to the CP and the voter model in the appropriate 
limits. The phase diagram of the model is obtained in 
one- and two-site approximations in section Hill In sec- 
tion IIVI results of simulations are shown which lead to 
numerical estimates of the critical fine and of dynamic 
critical exponents. Other estimates of the critical fine 
and static critical exponents are obtained through the 
exact diagonalization of the time evolution operator of 
the model for finite lattices, extrapolated to the thermo- 
dynamic limit using finite size scaHng in section^ Final 
comments and the conclusion may be found in section 
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II. DEFINITION OF THE MODEL 

Each one of the N sites of a one-dimensional lattice 
with periodic boundary conditions are occupied by par- 
ticles A or B. No holes are allowed. Thus, the state of the 
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system at a given time t is described by the occupation 
variables t] = ?72, Vn), where ry^ = or 1 if site i is 
occupied by particles B or A, respectively. 

The time evolution of the system is defined by the fol- 
lowing Markovian rules: 

1. A site i of the lattice is chosen randomly. 

2. If the site is occupied by a particle B, it becomes 
occupied with a particle A with a transition rate 
equal to Pana/2, where Ua is the number of A par- 
ticles in first neighbor sites of i. 
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Figure 1: Autocatalytic creation of a particle A at site i. 
Full circles represent A particles and empty circles denote B 
particles. 



3. If site i is occupied by a particle A, it may become 
a particle B through two processes 

• Spontaneously, with a transition rate Pc 



Figure 2: Example of spontaneous creation of particle B. 



Through an autocatalytic reaction, with a rate 
Pfcnh/2, where rib is the number of B particles 
in the first neighbors of site i. 
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Figure 3: Example of autocatalytic creation of a particle B. 

We define the time step in such a way that the param- 
eters Pa, Pb, and Pc, which are non-negative, obey the 
normalization Pa + Pb + Pc = 1, so that only two of them 
are independent. For convenience, we will discuss the be- 
havior of the model in the {pa,Pc) plane without loss of 
generality. 

The probability P{rj, t) to find the system in state ry at 
time t obeys the master equation 



dP{il,t) 
dt 



where r/* corresponds to the following configuration 

77* = (?7i,...,l -?7,,...,77Ar) (2) 



and Wi{ri) is the transition rate of the model, given by 



(3) 



where A =pa/{l - Pa), 1 = {I -Pc)/Pa, and the sum is 
over first neighbors of site i. 

The equation for the time evolution of the mean num- 
ber of A particles at site i, (rji), may be obtained from 
equations 1^ and being given by 



d{m) _ A 



dt 



(4) 



and a homogeneous solution, that is, p = (rji) and 
{rjirji+s) Vi, is given by equation 

dp 



^ = A[p - (2 - 7)</.] 



(5) 



In principle, we are not able to solve equation Q be- 
cause the function (j){t) is not known. We may write a 
differential equation for <^(i), but in this equation three 
variable terms such as {rjiijjrjk) will appear, so that a 
infinite hierarchy of equations will be obtained. A sys- 
tematic approximate solution is given below. Although it 
does not provide precise results, it still furnishes a qual- 
itative picture of the behavior of the system. 

This model, similarly to what happens in other out of 
equilibrium systems, displays absorbing states, which are 
such that once they are reached, they will never be left. 
The evolution rules of the model define the state where 
all sites are occupied by B particles (p = 0) as absorbing. 
Besides this stationary state, others may exist such that 
p — limt^oo pit) is nonzero. Such states are called active. 

It may be useful to remark that this model may be 
mapped to a spin system if we describe sites occupied by 
A and B particles by an Ising spin variables Ci = 1 and 
(Tj = —1, respectively. In these variables, the transition 
rate will be given by 



(6) 



where a = {pa + Pb + '2pc)/2, (3 = [pa ~ Pb - 'i.Pc)/ [pa + 
Pb + 2pc), e = {Pa + Pb)/{Pa +Pb + 2pc), and ^ = (p^ - 

Pb)/{Pa +Pb + '2pc)- 

Finally, this model corresponds to two known models 
in particular Hmits. If we make = or 7 = 1 the well 
known contact process is recovered 



w. 



(7) 



If now we take Pa = Pb and j3c = in the spin formulation 
of the model, the zero temperature linear Glauber model 
is recovered, also known as the voter model. 



w. 



(LGM) 



1 



(8) 
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III. MEAN FIELD SOLUTION 

An approximate solution of equation Q is obtained if 
we write (j) as a function of p. This approximation does 
not account for correlations and is known as the one site 
mean field approximation, assuming (?7i?7j) = 



that is 



may be written as 



Within this approximation, equation © 



| = (A-l)p-(2-7)Ap^ 



(9) 



and the stationary solution in terms of the parameters 
Pa and pc is 



P 







ifPa < 1/2 
otherwise. 



(10) 



Therefore, in the one-site approximation the absorbing 
state is separated from the active state by a continuous 
transition line located at Pa = 1/2 and the behavior of 
p as a function of Pa for some values of Pc is shown in 
figure 0. As expected, as pc is increased, the density 
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Figure 4: Density in the steady state as a function of pa for 
some fixed values of Pc . 

of A sites decreases, since annihilation of A particles is 
favored. It may be shown that close to the transition line 
the order parameter behaves as p ~ A'', A = {pa — 1/2) > 
0, with /? = 1 and, at the critical line, the density decays 
as p ~ t~'^, with ip = 1. These exponents are identical 
to the mean field results for the contact process. 

A somewhat better approximation is obtained if the 
set of equations for p and is simultaneously solved. 
As remarked above, an equation for (f) = (ijirjj) shows 
mean values of products of three rj variables. Therefore, 
in order to get a closed set of two equations, these mean 
values should be approximated as functions of (j) and p. In 
doing so, we are neglecting higher order correlations and 
this is known as the two-site mean field approximation. 

The following notation is helpful to obtain the second 
differential equation of the two-site approximation. We 
call p = where P{») is the probability that a site 



is occupied by an A particle and the probability to find 
a B particle is P{o) = 1 — P(»). We also define the 
three probabilities of the configurations of two neighbor- 
ing sites, and using the relations 



P(o) = F(.o) + P(oo), 
P(,) = P(.o) + P(..), 



(11) 

(12) 



together with the evolution rules stated above, the fol- 
lowing equation for P(»o) is found 

- \{Pa[Pioo.)-Pi.O.)] 



Pb[P{^ 



P(o.o)]} 



+ PcP{»»)-^{l+Pc)P{o»). (13) 

As already mentioned the probabilities of three site clus- 
ters appear in the equation. These will be written in 
terms of two-site probabilities through the so called pair 
approximation 



P{nin2nz) 



P(nin2)P(n2Ti3) 
P{n2) 



(14) 



Within this approximation, a closed system of two 
equations is found 



dp 
'dt 
du 
Hi 



= {pa-Pb)u-pcP 
,,2 



Pa- 



1 



Pb- 



SU - 



+ PcS - -(1 



(15) 



(16) 



where u = P(»o), r = P(oo) = I — p — u and s = P(»») = 

p — u. The first equation above corresponds to equation 
lO, whereas the second is obtained from equation lfT3jl 
using the pair approximation ifTHl . 

The stationary solutions for the densities are 

\{pa<p'a 

(17) 

Si-Zlt2;^P^+i-v. otherwise, 



where — |(4— Pc+ ^/SpTTpf). In this approximation 
the critical value of Pa is a function oi Pc- The critical 
exponents, as expected, have still their mean-field val- 
ues, only non-universal parameters are different in the 
two-site approximation as compared to the one-site cal- 
culation above. Figure © shows p as a function of pa 
for some values of Pc in the two-site approximation. 

It may be remarked that in the limit Pc ^ both 
approximations show a discontinuous transition between 
two absorbing states, which happens at pa = 1/2. In 
this limit the model corresponds to a spin model where 
a given spin is reversed with a probability proportional 
to the number of spins in first neighbor sites which are 
in the opposite direction. This might be called a biased 
voter model when Pa is different from pb . 
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Figure 5: Density as a function of Pa for some values of Pc 
results of the two-site approximation. 



IV. SIMULATIONS 

Although, as seen above, the model exhibits an ac- 
tive state in part of the parameter space in the approxi- 
mate solutions, numerical simulations are useful provid- 
ing some confidence that this state is not an artifact of the 
approximation. To our knowledge, for non-equilibrium 
models an argument similar to the Peierls construction 
for equilibrium models is still not available, and there- 
fore no systematic proof of the existence of the active 
state may be given, although this was accomplished us- 
ing probabilistic arguments by Harris in his original work 
on the contact process 

The simulation is done for discrete time, and may be 
described by the following steps: 

1. Initially, the system of TV sites with periodic bound- 
ary conditions is in a state with just one A particle. 

2. A list of all sites occupied by A particles is stored, 
and at each time step one of them is chosen ran- 
domly. 

3. Once the site is chosen, a random number p uni- 
formly distributed in the interval [0, 1] is generated, 
if p < Pa a B particle is replaced by an A particle 
in one of the first neighbors, if possible. Other- 
wise, the A particle at the chosen site will be turned 
into a B particle either through the spontaneous or 
through the autocatalytic process. 

4. To define the process switching A to B, another 
random number q is generated. If g < Pc/ {pb + Pc) 
the change is spontaneous, otherwise it will happen 
with a probability proportional to the number of B 
particles in first neighbors of the chosen site. 

5. The time interval associated with the steps above 
is At = 1/Na, where Na is the number of the sites 



occupied by A particles before the step. The pro- 
cess is repeated until either a maximum time t^ax is 
attained or the absorbing state Na = is reached. 

6. Several runs are done and mean values are calcu- 
lated as a function of time. 

The mean number of A particles, (Na), as a func- 
tion of time for simulations with — 10000 sites, 
tmax = 100000 and 10000 repetitions is displayed in fig- 
ure ^ . The number of sites is sufficiently high to ensure 
that in the simulations the cluster of A particles is much 
smaller than the system, thus avoiding finite-size effects 
in the time interval considered. In a region of the {pa,Pc) 
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Figure 6: Results of simulations 

plane an active state is found for long times in the simu- 
lations. Three curves are shown in figure ©. The convex 
one corresponds to an active stationary state, while the 
concave curve signals that the system will reach the ab- 
sorbing state. These two curves are separated by the 
third, where a power law (Na) ~ is seen, correspond- 
ing to the critical condition. We adopt, for all critical 
exponents, the notation proposed in the review article 
by Hinrichsen @. 

Repeating simulations in the region {pa + Pc < 1), 
the critical line may be estimated. The results of these 
calculations are shown in figure Q. The estimate for 
the critical line provided by the simulations display a 
concavity which agrees with the results of the the two- 
site approximation and, as expected, is situated above 
the mean field results, which usually overestimate the 
active region in the parameter space. The value of 
A = Pa/(1 — Pa) at the critical line attains the limiting 
value Ac = 3.2945 ± 0.0116 as pb 0, in good agreement 
with the estimated critical value for the one-dimensional 
contact process Ac = 3.29785(2) |il3|, as may be appreci- 
ated in figure (jHJ. 

The simulational result for (Na) as a function of time 
in the limit of the voter model {pa = l/2,pc = 0) is 
displayed in figure (EJ . One notices that the initial num- 
ber of particles is almost conserved, showing a narrow 
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Figure 7: Phase diagram of the model, showing estimates for 
the transition line between absorbing and active states in the 
one (MFl) and two (MF2) sites mean field approximations 
and simulations (S). In the line Pa + Pc = 1 the model cor- 
responds to the contact process (CP) and the unbiased voter 
model fLGMl is recovered at the Doint (vn = 1/2. Vr = 0) 
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Figure 8: Critical values of A = Pa/{1 — Pa) as functions of 
Pb- The dashed curve is a fit of a fifth degree polynomial to 
the points obtained in the simulations. 



dispersion. This aspect may also be appreciated in the 
histograms shown in figure (jlOII . for (pa = ^/2,Pc = 0) 
and {pa = 0.49, Pc = 0.001), with initially one A particle. 
For the voter model a small dispersion is shown around 
the maximum value (Na) = 1, while in the other case 
the maximum is shifted to the absorbing state {Na) = 

The simulations allow us to follow the dynamic evo- 
lution of mean values that describe the model, such as 
the mean number of A particles {NA{t)), the probability 
of survival Ps (t) that at the time t at least one particle 
A is present in the system and the mean square radius 
R'^it) = {T,z^'^Vi{t))/{NA{t)). These variables satisfy the 



Figure 9: Simulational results for the mean number of A par- 
ticles as a function of time for the voter model. 




Figure 10: Histograms of the mean number of A particles. 
The upper figure corresponds to the voter model {pa = 
1/2, Pc = 0) and the other figure is for a point of the param- 
eter space close to the voter model {pa = 0.49, Pc = 0.001) 
where the absorbing state {Na) = is stable. 



following scaling forms at the critical point 0: 



{NA){t) 

Ps{t) 



2/z 



(18) 



Through scaling and hyperscaling relations, these critical 
exponents define all other exponents of the model, and 
thus their values identify the universality class. 

From the simulations, the dynamic critical exponents 
defined above were estimated as functions of the rate Pc, 
as is shown in figure lfTT|l . For most of the values for Pc, 
particularly those close to the CP {pc ~ 0.2327), the esti- 
mated exponents are close to the values known for the CP 
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Figure 11: Dynamical critical exponents as functions of the 
rate Pc. The dashed lines indicate the values of these expo- 
nents for the CP. 



Figure 13: Verification of the hyperscaling relation for the 
dynamical critical exponents. 



(6 = 0.313686, 6 = 0.159464, and 2/z = 1.265226 
The relative errors are displayed in figure ltT2|l . Also, the 
hyperscaling relation 2d/z = 46 + 29 is satisfied numer- 
ically with some imprecision, as may be seen in figure 
II13II . Some systematic error seems to be present in the 
estimates for the mean square radius, which propagates 
to the hyperscaling relation verification. As pc 0, how- 
ever, when the voter model is approached, a departure of 
all estimates from the CP values is apparent. 
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Figure 12: Relative errors of the estimates of dynamical crit- 
ical exponents compared to the CP values. 



For models with a discontinuous transition with a spin 
inversion symmetry, such as the voter model, the dy- 
namic critical exponents are given by 0s = 0,5s = 1/2 
e 2/zs = 1 and the hyperscaling relation is given by 
6s + 9s = d/zs, corresponding to the compact directed 
percolation (CDP) universality class An analysis of 
the exponents in the limit Pc ^ may be performed fit- 
ting a cubic curve to the estimates, as is shown in figure 
itTljl for the exponent 6. The values obtained through 
this procedure are 6 = 0.44, 9 = 0.07, and 2/z = 1.01, 
which are not far from the known values given above. 
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Figure 14: Limit for pc = of the dynamical exponent 5. 



The variation of the estimates for critical exponents with 
the parameters of the model close to Pc = is proba- 
bly an apparent effect, due to large fiuctuations observed 
in the simulations in this region, and the whole critical 
line, with exception of the terminal voter model point, 
is in the same universality class of the CP, as suggested 
by the DP conjecture. This situation is quite similar to 
what is observed in the behavior of the Domany-Kinzel 
probabilistic cellular automaton (DKPCA) which in 
a line of its phase diagram corresponds to the CP with 
discrete time and parallel update, and also exhibits a 
crossover from the DP to the compact direct percolation 
(CDP) universality class It is believed, although 

to our knowledge not proven, that changing the update 
process in a stochastic model may affect non-universal 
quantities only. The point Pa = 1/2, Pc = may be rec- 
ognized as a multicritical point, and in its neighborhood 
any stationary density variable should exhibit the scaling 
form 



giPa-l/2,p,)^ipa~l/2r^F 



Pc 



(19) 



,ba-l/2]* 

The critical exponent associated with the density vari- 



able g, eg, should correspond to the CDP universality 
class, and the scaHng function F{z) is singular at a value 
zq of its argument, which corresponds to the critical 
line. Thus, the critical line is asymptotically given by 
Pc — zo{pa — l/2)'^ and the amplitude zq and the crossover 
exponent (j) may be estimated through a fit to the simu- 
lational results. Due to the large fluctuations mentioned 
above, these estimates are not very precise. The values 
we obtained are zq = 0.47 ± 0.08 and = 1.80 ± 0.03. 
In the two site approximation, we found <f> = 2 both 
for the model described here and for the DKPCA. Al- 
though extensive simulational results for the DKPCA are 
reported in the literature , apparently no estimate for 
the crossover exponent was obtained from them. 



V. NUMERICAL DIAGONALIZATION 
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Figure 15: Gap F as a function of Pa for system sizes A'^ = 
3,4,5,6 with pc = 0.2. 



Another approach to study 1+1 dimensional stochas- 
tic systems is the exact solution of models with increas- 
ing numbers of sites N followed by extrapolations to the 
thermodynamic limit This is accomplished writing 
the master equation as 



d_ 

dt 



P(t))^S\P{t)), 



(20) 



where S is the time evolution operator of the model. 
In the representation in which it is diagonal, vanishing 
eigenvalues = correspond to stationary states of the 
system. For finite systems with absorbing states, only 
these states are stationary, and no active stationary state 
is found. 

To study the transition between an active and a sta- 
tionary state, we may consider the behavior of the eigen- 
value with the second smallest absolute value F = /xi of 
the operator S. This eigenvalue is related to the quasi- 
stationary state and will eventually become degener- 
ate with fXQ in the thermodynamic limit, originating the 
phase transition. Figure II15|I shows the behavior of the 
gap F for some values of the size N . For a fixed value of 



which even for the limited number of estimates consid- 
ered lead to rather precise results. Figure lll(i|l shows 
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Figure 16: The squares are the extrapolated estimates ob- 
tained from the numerical diagonalization results, compared 
with the critical lines resulting from simulations (full line) and 
the two-site approximation (broken line). 



Pc- the finite size scaling behavior of the gap F is given 
by i 



T^N-^f{\pa-pl\N^/''-), 



(21) 



where is the critical value of Pa for a fixed value of Pc- 
Defining the quantity 



YN(pa,Pc) 



\n[T{pa,pc;N+l)/T{pa,pc;N~l)] 
ln[(iV+ 1)/(7V- 1)] 



,(22) 



we may estimate the critical point y'j {N) finding the in- 
tersection of the curves Yjsi and Yat+i [l^. This pro- 
cedure resembles the phenomenological renormalization 
group. The sequence Pa{N) of estimates for a given value 
of Pc, with = 4, 5, . . . M, was extrapolated to the ther- 
modynamic limit N ^ CO using the BST algorithm [lfl| . 



the extrapolated estimates for the critical line, compared 
with the results provided by the simulations and the two- 
site approximation results. The agreement between the 
simulation and the present results is apparent. The re- 
sults for the critical line close to the CP point ph — 
were extrapolated to the CP limit using a cubic function 
and the resulting value was Ac = 3.3081 ± 0.0173, which 
is in reasonable agreement with a more precise estimate 
in the literature Ac = 3.29785(2) We also estimate 
the crossover exponent <j), through a fit to the critical 
line in the neighborhood of the voter model limit and it 
is results in = 2.24 ± 0.07. 

Once the extrapolated estimates for the critical line 
were obtained, the critical exponents z and ^ = z — l/v±_ 
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may be estimated using the asymptotic scaling forms 

^T{pa,Pc) - L-i, (23) 

for a fixed value of Pc and Pa = Pa, on the critical line. 
An estimate of the exponents is obtained for each size 
of the system, and the estimates are extrapolated to the 
thermodynamic limit. Our results for these two static 
exponents are shown in figure l|17ll . 
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Figure 17: Static exponents estimated through the exact di- 
agonalization for finite systems. The values for the CP are 
shown in dashed lines. 

Again a departure of the estimates from the CP values 
may be observed as Pc becomes smaller. The CDP values 
for these exponents are Zs = 2 and = 1- Unlike to what 
was observed in the simulations no systematic trend to 
these values was observed in the estimates from exact 
diagonalization as Pc 0. We noticed, even considering 
values of N up to 15, that as Pc becomes smaller the 
convergence properties are poorer, and larger sizes may 
be necessary to obtain better estimates in this region. 
The estimates which were obtained from the simulational 
data seem to be more reliable that the ones which follow 
from the numerical diagonaHzation results. 



and ends at a point where the model is equivalent to 
the voter model, showing a discontinuous transition be- 
tween two absorbing states. The phase diagram thus 
obtained is similar to what is found in a problem of equi- 
librium polymerization in a grand-canonical ensemble on 
an anisotropic square lattice [20], where the polymer is 
modeled as a self-avoiding walk and the activity of a hor- 
izontal link is equal to x and a vertical link is associated 
with an activity y. A continuous transition is observed 
in general in the {x, y) plane, between a non-polymerized 
(p = 0) and a polymerized {p > 0), where p corresponds 
to the fraction of lattice sites visited by the walk. How- 
ever, at a; = or y = the walk is one-dimensional, 
and a discontinuous transition is found between the non- 
polymerized and the fully polymerized (p = 1) phases 

The localization of the critical line in the parameter 
space of the model was found in the one- and two-site 
approximations, as well as through simulations and nu- 
merical diagonaHzation of the time evolution operator. 
The evidences obtained from estimates of critical expo- 
nents suggest that the whole critical line belongs to the 
DP universality class, and a crossover to the CDP uni- 
versality class is observed at the terminal voter model 
point. From the estimates for the location of the critical 
line close to the voter model point the crossover exponent 
(f) — 1.80 ± 0.03 was obtained. Although the estimated 
value for (j) is smaller than the mean-field result (f) — 2, 
the latter is inside the confidence interval of the estimate, 
so that it is not possible to decide if the crossover expo- 
nent has a non-classical value. Due to the difficulties 
we had with the simulations and the numerical diagonal- 
ization in the multicritical region, we suspect that an- 
other approach is necessary to clear this point. We are 
presently addressing this point through a series expan- 
sion approach. Another extension of this work would be 
the study of the crossover between the universality class 
in the Domany-Kinzel cellular automaton, which corre- 
sponds to the CP with parallel update in a subspace of its 
parameter space and also has a point of its phase diagram 
which corresponds to the voter model. 



VI. CONLUSION 

We discussed the transition between an absorbing {p = 
0) and an active {p > 0) steady states in a stochastic 
model of interacting particles. This model displays a 
line of continuous transition between these states start- 
ing at a point which corresponds to the contact process 
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